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Growth cones of elongating neurites exert force against the external environment, but little is known about 
the role of force in outgrowth or its relationship to the mechanical organization of neurons. We used 
traction force microscopy to examine patterns of force in growth cones of regenerating Aplysia bag cell 
neurons. We find that traction is highest in the peripheral actin-rich domain and internal stress reaches a 
plateau near the transition between peripheral and central microtubule-rich domains. Integrating stress 
over the area of the growth cone reveals that total scalar force increases with area but net tension on the 
neurite does not. Tensions fall within a limited range while a substantial fraction of the total force can be 
balanced locally within the growth cone. Although traction continuously redistributes during extension and 
retraction of the peripheral domain, tension is stable over time, suggesting that tension is a tightly regulated 
property of the neurite independent of growth cone dynamics. We observe that redistribution of traction in 
the peripheral domain can reorient the end of the neurite shaft. This suggests a role for off-axis force in 
growth cone turning and neuronal guidance. 

Development of the nervous system is a complex series of events that involves large scale tissue morpho- 
genesis and migrations of cells. These events involve mechanical forces, but the mechanical properties of 
the nervous system are poorly understood compared to its chemical and electrophysiological properties. 
An exemplary feature of the nervous system is the great distance that neurites must extend from the cell body to 
their synaptic target'. At once each neuron must function as a single connected unit and as a set of compartments 
capable of responding to local stimuli. Abundant evidence indicates that neurites are under tension both in 
vitro^'^ and in vfvo" suggesting one means of cell-scale mechanical integration. The physiological significance 
of this tension is not fuUy understood, but several roles have been suggested". Tension can be a stimulus for 
growth during growth cone-driven neurite elongation', neurite towing by moving target cells'^, and the post- 
synaptogenesis stretch that elongates neurites during organismal growth' ''' '\ Tension can also specify axonal 
identity"", drive neurite pruning'^ and may contribute to compact wiring of the nervous system"* '''. 

It follows by simple force balance that growth cones, the dynamic, migrating tips of elongating neurites, must 
exert traction equal and opposite to neurite tension to resist retraction. Growth cone tractions have been 
measured using traction force microscopy (TFM) on compliant substrates in several recent studies^""^', however 
the mechanism of force generation is mostly unknown, as is its role in powering outgrowth. Traction force has 
been explored extensively in other motile cells types^*"^**, but growth cones present several unique features and 
experimental challenges. Growth cones exert much lower forces than other cell types consistent with their 
adaption to the relatively soft environment of the nervous system. Growth cones are also tethered to the cell 
body and other growth cones through the neurite shaft and therefore must be considered as a compartment of a 
much larger structure. The mechanical coupling between the growth cone and the neurite thus becomes crucial 
for understanding outgrowth. 

The growth cone is organized into compartments with distinct cytoskeletal architechture. The peripheral 
domain (labeled P, Figure la) is composed of a highly dynamic network of crosslinked actin filaments^"'. 
The central domain (labeled C, Figure la) contains a parallel array of uncrosslinked microtubules and their 
organelle cargos as well as stable bundles of actin filaments'"''". In the transition zone (labeled T, Figure la), the 
peripheral actin network condenses into contractile arc structures'^''', which are in turn connected to central 
domain actin bundles and the contractile actin node in the growth cone neck (labeled N, Figure la)" ''*. 
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Figure 1 | Spatial distribution of force in growth cones. ( a) Representative phase contrast image of Aplysia bag cell growth cone. P = peripheral domain, 
T = transition zone, C = central domain, N = growth cone neck. Scale bar is 20 /im. (b) Map of traction stress vectors in growth cone from (a) showing 
localization of traction force in peripheral domain. White outline shows edge of growth cone and neurite. (c) Map of strain energy density for growth cone 
in (a). Green line shows location of linescan. (d) Linescans of strain energy density through 38 growth cones. Traces are shifted such that the boundaries 
between peripheral and central domains align at zero. For each trace, values are normalized to the maximum value in that trace, (e) Internal stress 
measured along the same Kne scans as (d) normalized to maximum value in each trace. 



Aplysia californica bag ceil neurons have been used extensively to 
study the cytoskeletal dynamics underlying growth cone behavior on 
account of their large size and spread morphology. Here we use TFM 
with soft silicone substrates to explore the spatial and temporal pat- 
tern of force during bag cell neurite elongation. We show that neurite 
shafts are under tension and that tension is balanced by traction force 
in the growth cone peripheral domain. Growth cones can also exert 
forces that are balanced in other parts of the growth cone, which 
account for a variable proportion of total force. Although peripheral 
force patterns are extremely dynamic the tension is consistent over 
time. The progression of traction patterns during outgrowth suggests 
that a constant load continuously redistributes as adhesions are bro- 
ken and formed. These observations suggest a physical mechanism 
for neurite turning wherein assymetry of peripheral traction force 
reorients the tension and bends the end of the neurite toward regions 
of higher force. 

Results and Discussion 

Growth cones generate traction force in the peripheral domain. 

The mechanical engagement between growth cones and the extra- 
cellular environment drives the elongation of developing and 
regenerating neurites. Force measurement using compliant sub- 
strates was undertaken to extend our knowledge of the role of 
traction in growth. We have used displacements of tracer beads in 



a soft silicone gel to calculate the stress field, a vector map of force per 
unit area, and strain energy density, a scalar measure of the energy 
expended by the growth cone to deform the substrate. Figure lb 
shows a typical stress distribution and Figure Ic the corresponding 
strain energy density. 

Traction stresses applied by the growth cone on the substrate are 
centripetal, pointed away from the leading edge and toward the 
central domain. This follows the direction of flow of the underlying 
actin cytoskeleton. Actin filaments assemble at the leading edge and 
disassemble throughout the peripheral domaim'^"". The push of 
leading edge polymerization and pull of myosin motors within the 
network drives steady centripetal movement called retrograde 
flow'^'"'^". We found that highest stress and strain energy are 
consistently localized in the peripheral domain (Fig. 1 b,c, 
Supplementary Fig. SI). The spatial distribution of traction force 
was examined in 38 growth cones in which the boundary between 
peripheral and central domains was clear (dashed yellow line, 
Fig. la). Linescans from the leading edge, through the point of high- 
est force, and into the growth cone neck (green line in Fig. Ic) allowed 
comparison of stress patterns among a large number of growth 
cones. The linescans confirm that traction is primarily within the 
peripheral domain, however, the position of the peak of strain energy 
relative to the leading edge varies greatly between growth cones as 
does the width of the peripheral domain (Fig. lb,c,d. Supplementary 
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Fig. SI). Traction against the substrate must be balanced by internal 
stress on the cytoskeleton and other intracellular structures. Internal 
stress can be calculated as the integral of the traction stresses between 
each point on the line and the leading edge^'^ '*". We find that it rises 
throughout the peripheral domain, reaches a plateau in the transition 
zone, and continues at this level into the central domain (Fig. le. 
Supplementary Fig. SI). 

Consistent with the localization of traction force to the peripheral 
domain, cytochalasin B, which caps barbed end polymerization sites 
at the leading edge'^, led to collapse of the peripheral domain and 
cessation of traction force (Supplementary Fig. S2). The peripheral 
domain coUapse was reversible with drug washout. Jasplakinolide, 
which inhibits actin network turnover"", also led to peripheral retrac- 
tion and inhibited traction force (Supplementary Fig. S2). These 
results show that continuous assembly and disassembly of the actin 
network is necessary to maintain peripheral domain traction. 

Maximum traction stresses were typically a few Pascal (mean 90th 
percentile stress 6.8 +/— 2.8 Pa), confirming that neurons are 
among the weakest cell types'" We measure maximum internal 
stresses in the range of a few kilopascal (Supplementary Fig. SI) 
showing that although growth cone tractions are relatively weak 
the accumulated levels of cytoskeletal stress might be considerable. 
These results are generically consistent, but differ in detail, with 
previous reports of traction force in growth cones. The peripheral 
domain actin network includes parallel bundles called filopodia and 
branched meshworks between filopodia called veils'""''''. We clearly 
see the highest force under regions with veils and relatively low force 
under long filopodia (Fig. 1, Supplementary Fig. SI), contrasting with 
previous reports that have focused on force generated by individual 
filopodia extending beyond the veiP"'"'''. While others have shown 
highest stress at the leading edge^"'^', we find variable positioning 
of peak stress within the peripheral domain (Fig. lb,c and 
Supplementary Fig. SI). The pattern of internal stress reported here 
contrasts with published maps showing converging centripetal and 
outward directed forces near the transition zone^". We have found 
that internal force reaches a plateau in the transition zone and does 
not decrease into the growth cone neck, implying that this force is 
either balanced in other parts of the growth cone or transmitted to 
the rest of the cell as tension in the neurite shaft. 

Neurite tension is independent of area. Since mechanical equli- 
brium requires that the vector sum of traction stresses over the 
entire cell must equal zero, if the vector sum over the area of an 
individual growth cone is non-zero then this net force must be 
transmitted to the rest of the cell. We found that vector sums of 
traction stresses in growth cones give a net force in the direction of 
the neurite shaft (large green arrow. Fig. lb), which is equivalent to 
the tension transmitted to the rest of the celF"'^'. The scalar sum of 
stress over the growth cone area gives the total force that the growth 
cone is capable of generating. This commonly computed quantity"'* ''"' 
does not have a physical meaning, but when compared to the tension 
gives insight into the geometry of stresses in the cell. Mean tension 
was 3.1 +/— 1.6 nN, while the mean total force was 6.5 +/— 3.8 nN 
(n = 117) (Fig. 2a). The higher variability of total force can be 
accounted for by the observation that total force scales with 
growth cone area (r^ = 0.55, blue symbols) while the tension is 
very weakly correlated with area (r^ = 0.14, red symbols) (Fig. 2c). 
The increase of total force with size cannot be accounted for by 
higher stress levels in larger growth cones since the 90th percentile 
of the stress shows no correlation with size (Fig. 2b). The values of 
total force are all greater than the values expected if the area 
contained only background noise (Fig. Ic, green line). The values 
for tension reported here are consistent with reported values 
measured both with TFM^"'^' and other methods'' ''''°. 

The discrepancy between tension and total force can be under- 
stood as a consequence of the geometry of stress vectors in the growth 



cone. The ratio of tension to total force (referred to hereafter as the 
'tension ratio') gives the proportion of the force that is not balanced 
locally within the growth cone and therefore must be balanced else- 
where in the cell. This quantity is functionally equivalent to the 
'unbalanced traction' introduced to characterize force at ceU-ceU 
junctions''''''^. For an isolated cell or compartment, mechanical equi- 
librium requires zero net force and therefore a tension ratio of zero, 
whereas if all stress in the growth cone contributed to resisting neur- 
ite tension the ratio would be close to one (Fig. 2d). We found that 
growth cones with low tension ratios often have stresses near the 
neck partially balancing the rearward-directed stresses in more distal 
regions (Fig. 2f). While contractile stress associated with the cell 
margins appears to be a universal feature of adherent cells, such 
partially isolated contractile compartments have, to our knowledge, 
never been reported and might be a unique feature of neurons. 

Tension ratios for growth cones with area lower than the median 
(mean 0.6 +/— 0.2, n = 58) are significantly higher than tension 
ratios for growth cones with area higher than the median (mean 0.5 
+ /— 0.2, n = 57, P = 2e-5) (Fig. 2e). Thus, larger growth cones 
balance a higher proportion of force locally. 

Several possible physiological functions for neurite tension have 
been suggested'"'" ". But what is the function of the additional force 
that does not contribute directly to resisting tension? Investigations 
of neuronal development in vivo have shown that growth cones 
within nerve tracts are thin and streamlined but pause in their migra- 
tion and become complex and spread out when they reach 'decision 
regions''"'. This could be interpreted as a mechanism for exploring 
space at a point where positional information has to be interpreted, 
guidance cues sensed, and directional changes executed. The large 
spread growth cones in the current data set show locally balanced 
forces that are not directly resisting tension and which might be 
involved in mechanosensing, exploring space for guidance cues, or 
initiating a change in growth trajectory. 

Tension is stable during outgrowth. To investigate possible correla- 
tions between traction and outgrowth we sampled traction fields 
every 10 or 30 seconds in actively migrating growth cones. This 
timescale allowed us to resolve the transitions from one morpho- 
logy and stress pattern to another, but also allowed us to track growth 
over tens of microns. Accompanying movies of strain energy density 
overlain on DIG images show that the growth cone and associated 
traction patterns are extremely dynamic during outgrowth 
(Fig. 3a,b). In contrast, total force and tension, quantities that are 
integrated over the entire growth cone, are surprisingly consistent 
throughout the period of observation (Fig. 3c,d). This observation is 
reminiscent of the pattern observed for DRG neurons'". For the 
population considered the standard deviation of the tension was 
1.9 nN and the standard deviation of the total force was 5.3 nN. In 
comparison, the average standard deviation of the tension for 
individual growth cones over time was 0.8 nN, and the average 
standard deviation of the total force was 1.7 nN. Growth cones 
advanced steadily during observation and changes in growth rate 
were not associated with changes in force (Fig. 3c,d and Supplemen- 
tary Fig. S3). This observation is similar to the findings of Koch et 
al.^'. Further, average rates of advance do not depend strongly on 
tension, tension ratio, or area (Fig. 3e,f g). It should be noted that 
three growth cones with tension >4 nN also showed the slowest 
growth (Fig. 3e). The sample size is insufficient to establish 
whether there is a correlation, but this point deserves further 
investigation. 

The observed stability of tension over time suggests homeostatic 
maintenance within a narrow range. In fact, neurons have been 
shown to actively regulate tension to maintain a consistent value'"''''. 
The lack of correlation between tension and growth rate appears to 
contradict a body of evidence connecting externally applied tension 
with elongation'^'"". However, correlation between growth and 
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Figure 2 | Growth cone force scaling and local balance, (a) Histograms of tension and total force, (b) Plot of 90th percentile stress against growth 
cone area, (c) Plot of total force (blue symbols) and tension (red symbols) and expected total force value from noise only (green line) against growth cone 
area, (d) Example stress vector field of a growth cone illustrating a high tension ratio. Green arrow indicates direction of tension, (e) Histograms of 
the tension ratio for growth cones above (right) or below (left) the median area, (f) Example stress vector field of a growth cone illustrating a low tension 
ratio. Green arrow indicates direction of tension. For all plots N = 117 growth cones. Scale bars are 20 /im. 



internally generated tension, such as that imposed by growth cone 
tugging on the neurite, has never been clearly shown. Our data can be 
reconciled with the observations of tension-driven growth by a 
model in which the structure of the neurite continuously readjusts 
to maintain a set tension, elongating when the tension is greater than 
the set point and shortening when it is lower. Such a mechanism 
would require a tight feedback between growth cone advance, ten- 
sion generation, and structural remodeling. 

The independence of growth rate and growth cone size makes 
clear that the large, spread peripheral domain is not necessarily 
a 'paused' state, as suggested by in vivo studies'"". The lack of 
correlation between tension ratio and growth implies that stress 
field geometry does not limit growth rate. The fluctuations in 
tension (Fig. 2 c,d) appear smaller relative to the mean than those 
reported for DRG neurons, and the pattern contrasts sharply with the 
intermittent forces and movement observed in hippocampal neu- 
rons^'. Both studies however demonstrate that a highly dynamic 
and rapidly turning over structure can generate a consistent total 
output. 



Dynamic traction patterns balance a consistent tension. The steady 
tension must be maintained by a growth cone that is continually 
migrating and remodeling. Accompanying videos show that the 
force distribution changes rapidly as veils and filopodia protrude 
and retract (Supplementary Videos 1-4). How are the fluctuations 
in stress distribution coordinated to produce a consistent tension and 
total force? By examining the evolution of strain energy in different 
parts of a growth cone we find that traction dissipation in one area 
often corresponds to a rise in another, usually more distal, region 
(Fig. 4, Supplementary Fig. S4). Traction does not completely vanish 
in a region before starting to rise elsewhere, suggesting that a constant 
load' is being redistributed over the continuously remodeling growth 
cone structure. 

Traction is rarely distributed evenly throughout the periphery and 
is often concentrated into one or more 'hot spots'. We observe that as 
traction force shifts the central cytoplasmic domain bends towards 
areas of high force (Fig. 4a, Supplementary Videos 1-4). This is 
consistent with previous observations that a change in the direction 
of tension^ or traction on ligand-coated beads can reorient neurites'"'. 
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Figure 3 | Growth cone traction force is consistent during outgrowth, (a) DIG images and (b) strain energy density maps of a representative example 
growth cone advancing over 12 min. Scale bar is 20 ^m. (c) Growth cone position (red line, right axis), tension (blue line, left axis), and total force 
(green line, left axis) over time for the growth cone in (a), (d) Growth cone position, tension, and total force for another representative growth cone. 
Growth rate vs. (e) tension, (f) tension ratio, (g) growth cone area (N = 12 growth cones). 






Figure 4 | Traction rises and falls in different areas of the growth cone during outgrowth, (a) Overlay of DIG and strain energy density (green) at three 
timepoints. Solid lines outline the central domain at 0 (yellow), 240 (orange), and 570 s (red). Scale bar is 20 fim. (b) Montage of growth sequence 
corresponding to region 1 in (a), with 30 s interval. Scale bar is 10 ;Um. (c) Montage corresponding to region 2, 30 s interval. 
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The direction of tension (green arrow, Fig. 4a, Supplementary Videos 
1-4) is typically parallel to the direction of central domain pointing. 
The existence of neurite tension implies a mechanical connection 
between the peripheral domain and neurite shaft, and the association 
between peripheral traction, tension direction, and central move- 
ment (Fig. 4a) reveals tight mechanical coupling between of growth 
cone cytoplasmic domains during outgrowth and turning. This beha- 
vior of the central domain is reminiscent of the morphological 
changes that accompany the interaction oiAplysia growth cones with 
restrained ligand-coated beads or physiological substrates'""^'. 

The rapid fluctuation of stress distribution is consistent with the 
current understanding of growth cone adhesion dynamics. Growth 
cones form highly dynamic structures called point contacts with a 
similar composition to the more stable focal adhesions of fibroblasts 
and other cells^^'^l These contacts are likely responsible for the mech- 
anical coupling between cytoskeleton and substrate. Leading edge 
protrusions are stabilized by the formation of point contacts^^ and 
modulation of contact dynamics may underly turning in response to 
guidance cues^*'^"". Point contacts have an average lifetime of 2.2 
minutes"'"', consistent with the time scale of appearance and dissipa- 
tion of traction in different regions of the growth cone. According to 
a popular model of growth cone motility^'''", traction between mov- 
ing actin and substrate-anchored adhesions drives outgrowth by 
resisting retrograde flow and converting polymerization into leading 
edge protrusion. Examination of the accompanying movies shows 
that leading edge advance and appearance of traction are often out 
of phase (see Supplementary Video 2, 360-540 seconds and 
Supplementary Video 3, 390-600 seconds), suggesting that the forces 
measured here might not be directly involved in driving protrusion. 
Simultaneous visualization of adhesion, actin, and traction dynamics 
will clarify the connection between cell-substrate coupling and 
growth cone advance. 

We propose a simple mechanical model for growth cone advance 
wherein the tension of the neurite shaft presents a constant load that 
is balanced by growth cone traction force (Fig. 5). This force is dis- 
tributed among cell-substrate adhesions under the peripheral 
domain. When some adhesions break the force becomes redistribu- 
ted among the existing adhesions. In our hypothesized model, this 
movement under constant load can be likened to a rock climber 
whose constant weight under gravity is continuously redistributed 
among shifting handholds and footholds. Thus growth may proceed 
by the formation of adhesions under the extending peripheral 
domain with subsequent loss of adhesions near the central domain 
and at the sides and rear of the growth cone. Our model also suggests 



a mechanism by which off-axis traction force could initiate neurite 
turning. In response to guidance cues or substrate borders, the select- 
ive stabilization of adhesions on one side of the growth cone could 
lead to a reorganization of traction force, a shift in the direction of 
tension, and turning of the neurite toward the attractive cue (Fig. 5 a 
to b and c to d). 

In conclusion, our measurements of growth cone traction force 
during neurite elongation have revealed that traction in the peri- 
pheral domain balances tension in the neurite shaft. This tension is 
independent of growth cone size while the total force increases with 
area. These data indicate that growth cones are mechanically inte- 
grated with the rest of the cell through neurite tension but can also be 
partly mechanically isolated when a portion of the force is balanced 
locally. The existence of partially isolated compartments may be a 
unique feature of neurons, necessary on account their extended mor- 
phology in which different parts of the cell must respond indepen- 
dently to their local environment. We also show that tension and 
total force, integral quantities describing growth cone force output, 
are consistent over time despite highly dynamic peripheral domain 
stress distributions. 

The narrow range of neurite tensions within the population and 
the consistency of tension over time suggests that this might be a 
tightly regulated property. Identifying the feedback mechanisms and 
signaling pathways that contribute to force regulation will be an 
important subject of future study. Limits to tension may not arise 
from the growth cone itself since large growth cones are clearly 
capable of generating tensions much larger than those measured, 
but might be a property of the neurite shaft. We propose that the 
growth cone moves under a constant load imposed by the neurite 
shaft. This load is continually redistributed among adhesions that 
form and break during advance. 

Methods 

Traction force substrate preparation. Traction force substrates were prepared by 
coating a glass coverslip with a highly elastic silicone gel as described^'*, with the 
following modifications. One layer of fluorescent beads {100 nm or 40 nm yellow- 
green or red fluorescent carboxylate-modified microspheres, Invitrogen) at a 1 : 2000 
volume ratio in 100% ethanol were deposited on the surface of an acid- washed 
coverslip by spin-coating at 2,000 RPM for 30 sec with a PWM32 spinner {Headway 
Research, Garland, TX). Silicone gel was prepared as a 1 : 1 weight ratio of CY52-276A 
and CY52-276B (Dow Corning Toray Silicones, Tokyo, Japan). The mixture was 
degassed for 10 minutes and applied to the coverslip by spin-coating for 30 seconds at 
2,000 RPM. The film of silicone gel prepared in this manner was typically between 30 
and 40 /tm thick. The silicone was cured overnight at room temperature. Silicone gels 
prepared in this manner had a Young's modulus of 3 kPa. 

After curing, fluorescent beads were deposited on the top surface of the silicone gel 
using a reaction of {3-Aminopropyl)triethoxysilane (APTES) and l-ethyl-3-{3- 



traction stress 




Figure 5 | Model for growth cone advance. Peripheral domain traction (red arrows) balances a constant load imposed by the neurite (orange arrow). 
When a portion of the growth cone retracts (a to b) the tension becomes redistributed among remaining points of attachment, accompanied by the 
movement of the central domain towards the region of higher force on the right. When a new region of peripheral domain advances (b to c) force is 
redistributed among the old and the newly formed adhesions. Detachment on right and advance of the peripheral domain on the left (c to d) again leads to 
the redistribution of traction and reorientation of the central domain. 
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dimethylaininopropyl)carbodiimide (EDC). APTES { Sigma- Aldrich, St. Louis, MO) 
was applied to the silicone gel surface by vapor deposition in a vacuum chamber for 4 
hours. 1 mg/mL EDC {Sigma- Aldrich) was dissovled in 50 mM borate buffer and 
100 nm green-yellow or red beads were added to a fmal dilution of 1 : 1,000 or 
1 : 2,000 and sonicated with a VirSonic 100 probe sonicator (VerTis Co., Gardiner, 
NY). Silicone gel coated coverslips were then floated on top of this solution for 10-15 
minutes. 

Cell culture and media. Prior to plating bag cell neurons, vertebrate laminin {Sigma- 
Aldrich) diluted with Aplysia hemolymph to 50 /(g/mL was passively adsorbed onto 
the surface of TFM substrates for 1 hr at room temperature. Hemolymph was 
collected as previously described^^. Substrates were rinsed with L15 artificial sea water 
(ASW) {400 mM NaCl, 10 mM KCl, 15 mM HEPES, pH 7.8, 10 mM CaCls, 55 mM 
MgCl2, and Phenol Red). Primary culture of Aplysia bag cell neurons was as 
previously described^^. Bag cells were grown overnight in L15-ASW either at room 
temperature or at 18 C prior to use in experiments. Experiments were performed at 
room temperature in L15-ASW supplmented with 0.5 mM vitamin E (Sigma- 
Aldrich) and 2 mg/ml carnosine (Sigma- Aldrich). Solutions containing Cytochalasin 
B {Sigma -Aldrich) and Jasplakinolide {Molecular Probes, Eugene, OR) were applied 
to cells using a homemade chamber in which a compartment containing the cells is 
separated from a solution exchange compartment by a Whatman Cyclopore track 
etched membrane with 5 /fm pores (Whatman, Florham Park, NJ). This arrangement 
was necessary to prevent mechanical disturbance of the cells during solution 
exchange. 

Image acquisition. Images were acquired on a Nikon TE2000E inverted microscope 
(Nikon, MelviUe, NY) and an Andor Revolution spinning disk confocal system 
(Andor, Belfast, UK) equiped with a CSU-Xl confocal head (Yokogawa, Tokyo, 
Japan) and a Andor iXon^'^"^888 EM CCD camera. Confocal illumination was with 
488 and 561 laser lines controlled with an Andor Laser Combiner. Emission 
wavelength was selected with bandpass filters from Chroma Technology (Bellows 
Falls, VT) mounted in a Sutter LBlOW-2800 filter wheel. Transillumination was with 
a halogen lamp and SmartShutter (Sutter Instruments, Novato, CA). AH hardware 
and image acquisition were controlled with /(-Manager software^". Objectives used 
were Plan Apo VC lOOX/1.4 numerical aperture (NA), phase contrast Plan Apo 
lOOX/1.45 NA, Plan Apo 60X/1.4 NA, and phase contrast Plan Apo 60X/1.4 NA 
(Nikon). 

TFM data consisted of single fields at single time points or a single field imaged at 
multiple time points. For each field a phase contrast or DIC image was acquired of the 
growth cone, then a confocal image was acquired of the upper layer of fluorescent 
beads. The growth cone was centered in the field of view such that bead displacements 
would be expected to decrease to zero at the edges. For all TFM experiments an 
unstressed reference image was required to calculate absolute stress magnitudes. 
Accordingly, after the experiment, chambers were perfused with 1% Triton X-100 
(Sigma- Aldrich) in ASW to dissolve cells and relax the stress exerted by the cells on 
the substrate. A confocal z-stack of the beads was then acquired for every field of 
interest. 

Calculation of traction fields, strain energy, and internal stress. AH TFM data was 
analyzed in MATLAB (MathWorks) and imaging processing was with either 
MATLAB or ImageJ^\ Bead positons were found using MATLAB's centroid 
analysis^^ from either single confocal images or average intensity projections of z- 
stacks. The deformation of the substrate, u,{r), was found by calculating the distance 
between centroids of stressed images and the unstressed reference image^^. Linear 
elasticity theory relates the traction stress to the deformation in Fourier space by (7,z(k, 
K) ~ Q(7{ki ^s)w,(k), where k is the in-plane wave vector and is the thickness of the 
substrate. Gi^O^, and t/y(k) are the Fourier transforms of the traction stress and bead 
displacements, respectively. The tensor Q is dependent on the elastic properties of the 
silicone gel, the wave vector, and /i/*'^^. The strain energy density, a scalar quantity 
equivalent to the energy per unit area expended by the cell to deform the substrate, 

was calculated as w(r) — -(Tiz{r,hs)ui{ry^. Tension was calculated as the magnitude 

of the vector sum of stresses over the area of the growth cone, | j dA (cr^z^ + ^^yzf) | > 

and the total force was calculated as the scalar sum of stress magnitudes over the area 

of the growth cone, j (M|(TjczX-|-(T^zy|> where x and y are unit vectors. Noise was 

estimated by finding the average stress magnitude, iiA|»T^2X-|-(7y2y|^ /A, for 

regions outside the growth cone in a subset of traction force maps. 

Linescans were made by drawing a line from the leading edge, through the region of 
highest strain energy, and into the neck of the growth cone. For strain energy density 
traces, the strain energy density map was sampled at regular intervals along the line. 
Stress along the line was found by projecting adjacent stess vectors onto the line. 
Internal stress along the line was estimated using o^x — I dxio^Jh) as described^^'*". 
The thickness of the growth cone, h, was taken to be 180 nm in the peripheral domain 
and 1300 nm in the central domain^^. The internal stress calculation assumes that all 
stress is oriented along the line and the gradient of the stress perpendicular to the line 
is small. This is a good approximation of our data since the linescan is parallel to the 
direction of stress and through a maximum. Growth was measured from DIC images 
as displacement of the growth cone's central domain along the inferred growth axis, 
sampled every 50 or 150 seconds. 



Statistical analysis. All statistical analyses were performed in MATLAB. p values 
were determined with a Wilcoxon paired-sample test for comparison of drug 
treatments to control, and with a two-sided Wicoxon rank sum test for comparison of 
tension ratios. Statistical significance was defined as p < 0.05. 
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